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Quantum interference of correlated particles is a fundamental quantum phenomenon which carries 
signatures of the statistics properties of the particles, such as bunching or anti-bunching. In presence 
of particular symmetries, interference effects take place with high visibility, one of the simplest cases 
being the suppression of coincident detection in the Hong-Ou-Mandel effect. Tichy et al. recently 
demonstrated a simple sufficient criterion for the suppression of output events in the more general 
case of Fourier multi-port beam splitters. Here we study the case in which 2 q particles (either 
bosonic or fermionic) are injected simultaneously in different ports of a Sylvester interferometer 
with 2 P > 2 q modes. In particular, we prove a necessary and sufficient criterion for a significant 
fraction of output states to be suppressed, for specific input configurations. This may find application 
in assessing the indistinguishability of multiple single photon sources and in the validation of boson 
sampling machines. 


I. INTRODUCTION 

Multi-particle quantum interference arises when sev¬ 
eral indistinguishable particles have non-vanishing prob¬ 
ability amplitude of being found at the same site or spa¬ 
tial coordinate. The algebraic sum of all these proba¬ 
bility amplitudes may lead to strong enhancement (con¬ 
structive interference) or suppression (destructive inter¬ 
ference) of the detection probability of the different pos¬ 
sible collective states. In that, it is a pure and typical 
quantum phenomenon, which is worth being investigated 
both from a fundamental perspective and for its quantum 
information implications. 

Qualitatively different behaviours may be observed in 
general, depending on the bosonic or fermionic nature of 
the particles. The anti-symmetrization requirements [I] 
for fermionic wavefunctions lead to vanishing probability 
of finding more than one particle on the same site (Pauli 
principle). Bosons, on the contrary, show a remarked 
tendency to bunch together, with increased probability 
of coalesce on the same site mi or to cluster in nearby 
sites (bosonic clouding [3]). 

However, when particles evolve following Hamiltonians 
with specific symmetries, particular fine-grained distribu¬ 
tions can be observed with enhanced interference peaks 
and dips. The simplest case is the Hong-Ou-Mandel 
(HOM) effect when two particles impinge on distinct 
ports of a balanced beam splitter: quantum interference 
suppresses the coincident output (one particle per each 
output) in the case of bosons and the single-port output 
(both particles in either output) in case of fermions. In 
the multi-particle case, a sort of generalised HOM effect 
occurs for symmetric multiport beam splitters [SHS]. In 
particular, Tichy et al. [7j showed that for a particular 
class of multiports, namely Bell or Fourier multiports, 
and input states with cyclic symmetry, a full suppres¬ 
sion of most of the output combinations is observed; a 
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simple analytical law gives a sufficient criterion for such 
suppression. 

From a computational point of view, calculating the 
output distribution of a number of indistinguishable 
bosons is a mathematically hard problem, in that it can¬ 
not be performed efficiently on conventional (classical) 
computers. In fact, it relies on the calculation of per¬ 
manents of matrices, for which an efficient classical algo¬ 
rithm is lacking. The realisation of such difficulty has led 
to the proposal [9] of boson sampling devices as exper¬ 
imentally accessible platforms that could perform some 
task hard-to-simulate with classical resources. The spe¬ 
cialised task of such quantum devices is to physically im¬ 
plement and sample the distribution of n bosons under¬ 
going a certain unitary evolution. The computational 
difficulty of a classical simulation of such process (i.e. 
a classical sampling of such distribution) increases expo¬ 
nentially with n , rapidly becoming infeasible. First proof 
of principle experiments with photons have been reported 
very recently mmm-, while not having demonstrated 
yet a true quantum supremacy, these experiments have 
pointed out that such a demonstration may not be so far. 

A future many-modes boson sampling experiment will 
likely require the implementation of an arbitrary unitary 
matrix, through a possibly reconfigurable m linear in¬ 
terferometer [SUE]. If the output is hard to predict 
classically, it may be not trivial also the certification of 
the correct operation of such device: in fact, several so¬ 
lutions to this problem have been debated HIHH7HE]. 
The use of particular symmetric unitaries that show rich 
but easily predictable multiphoton distributions has been 
also proposed [18] as a convenient way to assess both the 
performance of a multiphoton source, as well as of the 
reconfigurable device itself. In the same way the two- 
photon HOM effect on balanced beam-splitters has long 
been used as a diagnostic instrument for measuring two- 
photon indistinguishability, suppression laws for multi- 
port interferometers could provide a suitable means to 
simultaneously test the quality of a multi-photon source 
and of a multi-mode reconfigurable device for boson sam¬ 
pling experiments. Of course, these could be adopted to 
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assess the performance of a multi-photon source also out¬ 
side of the boson sampling context. 

It has to be noted that the existence of a sharp sup¬ 
pression law directly comes from the symmetry charac¬ 
teristics of the matrix. While it has been conjectured 
that other suppression laws could exist for other class of 
symmetric unitaries, only the class of Fourier matrices 
has been investigated extensively up to now[Tl BUI [21] , 

In this work we address the study of interferome¬ 
ters implementing m = 2 p -modes Sylvester matrices and 
prove a necessary and sufficient criterion for the suppres¬ 
sion of most output combinations, for certain input states 
of n = 2 q < m particles, either bosonic or fermionic. We 
further discuss possible applications for assessing the in- 
distinguishability of multiple single-photon sources. 

In Section [TT] we recall some basic concepts about the 
evolution of multi-particle Fock states through linear uni¬ 
tary processes and about the definition of Sylvester ma¬ 
trices. In Section III we give a comprehensive charac¬ 
terization of the output distributions in the two-particle 
case, while in Sections [TV] and [V] we prove the output 
suppression criteria for the cases of multiple bosons and 
multiple fermions respectively. These sections are mainly 
organised as a list of propositions regarding mathematical 
properties of certain matrices and their consequences on 
the calculation of the multi-particle output distributions. 
Finally, we discuss in Section |VI| the consequences and 
possible applications of the suppression criteria proved 
in the preceding sections, with particular regard to ex¬ 
periments with photons. 


II. PRELIMINARY CONCEPTS 


A. Multi-particle interference in linear 
interferometers 


A generic Fock state | T) of n particles on m modes can 
be written as | T) = ^Iir=i a t,) | 0 ) where aj. is the cre¬ 
ation operator on the mode t,. Such state can be iden¬ 
tified by the ro-element vector t = (ti, t 2 , ■ • •, t n ), with 
1 < U < m. Since different orderings of the particles in 
the same modes are not distinguishable, we will consider 
only the cases t\ < £2 < ... < t n . 

An ?n-mode lossless linear evolution can be described 
by a to x to unitary transformation U on the space 
of creation operators. The probability amplitude as¬ 
sociated with an input g = (gi, g 2 , ■ ■ ., g n ) and output 
h = ( hi , h 2 , ■ ■ ■, h n ) is given by 


= _ P eimS g,K _ 

Pb0b ■ ■■Hm'-VlW- ■ • ■ I'm! 

in the case of bosons and by 


Pfer — 


det S-. r 

_ 9jh _ 

V/rW ■ • ■ Vm)-V l!zx 2 ! . . . V m \ 


(1) 


( 2 ) 


where gi and v- L are the number of particles present in 
mode i in the g and h states respectively, S is the scat¬ 
tering matrix with elements Sij = Uh it3j and perm A 
denotes the permanent of a matrix A. It may be useful 
to recall also the definition of permanent for a matrix A: 

n 

perm A = ^ a i>er (i) (3) 

a i— 1 

where aij is an element of A, a is a permutation of 
{ 1 ,..., n} and thus the sum in the expression is per¬ 
formed over all the possible permutations. 


B. Hadamard, Fourier and Sylvester matrices 


A complex Hadamard matrix is defined as an orthogo¬ 
nal matrix of complex numbers, in which all the elements 
have unitary modulus. A well-known sub-class of such 
matrices is that of Fourier matrices, the elements of a 
to x to Fourier matrix F{m) being defined as follows: 

F jtk = e 2^0-i)(fc-i)/™ ( 4 ) 


where i — y/—T is the imaginary unit. As already men¬ 
tioned, multi-particle interference has been largely stud¬ 
ied in the literature [7J [201 [2T| for interferometers imple¬ 
menting the normalized (unitary) version of such matri¬ 
ces U m = -^=F(m). 

Real Hadamard matrices, simply referred to as 
Hadamard matrices in the following, are orthogonal ma¬ 
trices with all elements equal to ±1. Sylvester matrices 
are a particular class of real Hadamard matrices, hav¬ 
ing size to = 2 P , that can be built recursively from the 
following formula: 


H(2 p ) = 


H(2 p - 1 ) 


H{2 p - 1 ) 

-H(2P~ 1 ) 


(5) 


being H( 2°) = H(l) = [1]. From this construction one 
can derive an analytic expression for the (i,j) element of 
the matrix 


[H(2 p )\ ij = (-1 (6) 

where is and j b are the binary representations of i and 
j, enumerating the rows and columns starting from 0 , 
and 0 is the bitwise dot-product. 

In the following we will refer to devices implementing 
a unitary matrix of the kind: 

U m = -L H(m) (7) 

v TO 

with to = 2 P as Sylvester interferometers. 

A general expression for the permanent of Sylvester 
(and more in general, Hadamard) matrices is not known: 
while perm if ( 2 ) = 0 , it has been conjectured that for all 
the other orders Hadamard matrices have non-vanishing 
permanents [ 22 ] . 
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Note that the usual balanced beam-splitter operator is 
just ~^H( 2). In this case, for an input state with one 
photon per mode, the well known Hong-Ou-Mandel ef¬ 
fect is observed, which consists in the suppression of the 
output state with one photon per mode. In fact, accord¬ 
ing to 0 , for such output contribution, the probability 
amplitude is proportional to permU( 2 ) = 0 . 


III. TWO PARTICLES 


Even though some of the results of this section could 
be retrieved by applying the more general results of Sec¬ 
tions |IV| and |Vj the two-particle case allows for a more 
comprehensive description and shows some specific fea¬ 
ture, which make it worth addressing it separately. 


Proposition 1 If two bosons are injected in the first 
two modes of an interferometer described by U m = 
H(2 P ), the probability amplitude Pij of 
an output state with one particle on mode i and one par¬ 
ticle on mode j follows the rule: 


\Pid\ - 2 P-1 

= 2 p-i/2 

Pi,j = 0 


mod 2 = j mod 2 , i ^ j 

if i = j 
else 


The scattering matrices 5 in 0, for such an input 
state, are all sub-matrices of the first two columns of U m . 
Reminding ([5]), one can easily observe that such columns 
are just (properly normalized) repetitions of H( 2) and 
retrieve the matrices 5, as a function of i and j, as fol¬ 
lows: 


• i mod 2=1 and 

S = 73?-^(2) = 775J 


j mod 2 = 0 

1 1 ~ 


72 p 


1 -1 


perm 5 = 0 


• i mod 2 = 0 and 
'1 -1 


S' - 1 

* “ 72 ^ 


1 1 


j mod 2 = 1 
perm 5 = 0 


• i mod 2 = j mod 2 

5 consists of two identical rows with elements 
± 7 j=: simple calculations show that in this case 

| perm 5| = 


Thus, if and only if i mod 2 mod 2 the permanent 
of the scattering matrix vanishes, giving pij = 0 , while 
the other cases are proved by applying 0 . □ 


Corollary 1.1 In the case of Prop. [7J the fraction of 


N° 


suppressed states is N B f vp = 2 ++ 


The possible two-bosons output states are identified by 
all the couples (i,j) with 1 < i < j < m (we consider 
only i < j because (i, j) is the same state as ( j,i )). The 


number of such states is ^to(to + 1). In a chessboard 
with to x to squares, alternately black and white, these 
can be seen as all the squares above the main diagonal 
or included in it. If we colour the squares in such a way 
that the main diagonal is black, the condition i mod 2 ^ 
j mod 2 (with i < j) indicates all the white squares 
comprised in the region above it, which are actually half 
of the total number of white squares. Thus, the number 
of suppressed states is \rv ?, giving the result. □ 

Proposition 2 If two fermions are injected in the first 
two modes of an interferometer described by U m = 
H(m ) = ^=iJ( 2 P ), the probability amplitude ptj of 
an output contribution with one particle on mode i and 
one particle on mode j follow the rule: 

I Pi,j I = if * m od 2 ± j mod 2 

Pij = 0 else 


The scattering matrices to be considered for calcu¬ 
lating the probability amplitudes are just the same of 
Prop. [TJ but the determinant (Eq. 0) instead of the 
permanent has to be calculated here. Thus, when i 
mod 2 7 ^ j mod 2 the scattering matrix is composed of 
two identical rows and the determinant vanishes. In all 
the other cases (see the expression of 5 in the proof of 
Prop, [l]) the determinant is equal to g+r- Application 
of (0 then gives the probability amplitudes. □ 


Corollary 2.1 In the case of Prop. \^the fraction of sup- 

7vrfer 

pressed states is » f supp = 

^ N b tates 2 m + 1 


Comparing Prop. [T] with Prop. [2] bosons and fermions 
show a dichotomic behaviour, in that an output combina¬ 
tion is suppressed for two bosons if and only if is allowed 
for two fermions and vice versa. Thus the fraction of sup¬ 
pressed states (over all the possible two-particle states) 

for two fermions is = 1 - = 1- b = 

-‘’states -‘’states 2 771+1 

1^ 771+2 I—I 
2 m+1 ’ LJ 

Note that this fraction actually includes some states 
(the states with two particles on the same port) that are 
indeed suppressed by virtue of the Pauli principle and 
not by specific features of the Sylvester matrix. 


Proposition 3 For U^p = ^=If(2 P ) and a two-particle 
input on an arbitrary couple of different modes, the num¬ 
ber of suppressed states is the same as that given in 
Corol. 03 and Corol. \2.1\ for bosons and fermions re¬ 
spectively. 


For an arbitrary input state, with two particles on 
modes (*,i), the scattering matrices 5 will take elements 
from the *-th and j-tli columns of ~^==H{ 2 P ). Let us put 
these two columns one next to the other, to form the n x 2 
matrix A. For a given output state ( i',j'), the scattering 
matrix 5 will be a submatrix of A formed by its th and 
j'-th rows. 
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Note that two different columns of an Hadamard ma¬ 
trix have half of the elements with opposite sign and half 
of the elements with the same sign. Half of the rows of 
A will be [1,1] or [-1,-1]; the other half will be [—1,1] 
or [1,-1], Let’s now perform the following operations. 
First, we multiply the [—1, —1] and [—1,1] rows by -1. 
This will change the sign of the permanent of the scat¬ 
tering matrix that should include such rows, but it has 
no influence if the permanent vanishes. At this point we 
will have half of the rows equal to [1,1] and half equal 
to [1,-1]. Second, we reorder the rows alternating [1,1] 
to [1,-1]. This is equivalent to relabelling the outputs, 
which does not affect the number of suppressed output 
states. At this point the matrix A will be just the same as 
if (i, j) = (1,2), which is the case of Prop, [l] and Prop. [2] 
Hence, every input combination has the same fraction of 
suppressed output states of the input (i,j) = (1,2), dis¬ 
cussed in Corob fOl and Corol. lAll for boson and fermions 
respectively. □ 

IV. MULTIPLE BOSONS 

The aim of this section is to demonstrate a suppression 
law for the case of n = 2 q indistinguishable bosons. This 
will be obtained in Prop. [7] at the end of the section. 
However, that result is based on other propositions which 
will be proved before. The first one (Prop. [4]) allows to 
restrict the study, in certain conditions, from the case 
of n particles in m modes to the case of n particles in n 
modes. Props. urn on the contrary, regard mathematical 
properties (in particular, the value of the permanent) of 
certain —1,+1 matrices. 

Proposition 4 Be U = ~^H(2 P ) with p = k + q a 

linear transformation over m = 2 P modes, and h = 
(1 + n ■ c,..., n + n ■ c), where 0 < c < (2 fc — 1), an 
input state of n = 2 q particles. The output state g = 
(< 7 i,... ,g n ) is suppressed if and only if the output state 
9 1 = (g[, ■ ■ ■, g' n ) with g\ = [(gi - 1) mod n] + 1 is sup¬ 
pressed for the transformation U' = H(2 q ) with n 
particles entering one per each mode. 

For input states of the kind h = (1, ...,n) (i.e., one 
particle per each of the first n modes), the scattering 
matrices will be submatrices of the first n columns of U. 
From the construction of H(2 P ) = H(2 k+g ) = H(n ■ 2 k ) 
with |5]), it is clear that these first n columns are just 
repetitions of H(n). For an output state with n parti¬ 
cles distributed on the modes g\.. .g n , the i-th row of 
the scattering matrix S will be extracted from the gi~ 
th row of U. Since such rows repeat identically every n 
rows (for that regards the first n columns), the scattering 
matrix is the same for all states g = (g 1 ,... ,g n ) having 
the same g 1 mod n,... ,g n mod n. We can look for the 
smallest mode numbers giving this condition, which are 
9i = Kdi ~ 1) m °d n] + 1. In that case the scattering 
matrix is the one we would expect for the output state 


g' = (g[,... ,g' n ) defined as above, when entering with n 
particles (one per each mode) in an interferometer imple¬ 
menting U' = 2 g ). 

Let’s consider now the more general case h = 
(1 + n • c,..., n + n ■ c) with 0 < c < 2 k ~ 1 . Again from 
the construction in ([ 5 ]) it can be observed that such 
columns will be repetitions of ±H ( n ) (with a succession 
of signs + and - that depends on c). If we properly change 
the signs of the rows (operation that is equivalent to add 
a 7 r phase term to certain outputs, which does not influ¬ 
ence the probability modulus) these columns can be made 
identical to those of the case h = (1, ...,n), discussed 
above. Hence, the output distribution is the same □. 

This result hold for both bosons and fermions be¬ 
cause no hypotheses on the particles statistics have been 
adopted. In addition it can be exploited for a more pre¬ 
cise generalisation of the results of Props. [T] and [2] to a 
wider range of input states. 

Proposition 5 Be A an m x m matrix, with m = 
2 P , built by taking the rows {r\, r 2 ,..., r m } from H(m) 
(namely, the i-th row of A is the ri-th row of H(m) and 
rows may be repeated). If q ^ ® ® • • • ® 7 ^ 0, 

then perm A = 0, being r^^s the binary representation of 
the row number, starting the count from 0, and ® the 
bitwise sum (XOR operation). 

The condition rqs ® T 2 ,b ffi • ■ ■ ® r m ,B 7 ^ 0 means that 
for at least one k , an odd number of r^B has the same 
fc-th bit. In other words, for at least one k , the fc-th bit 
of the binary representations of the r,; is 1 for an odd 
number of rows and is 0 for a (possibly different) odd 
number of rows. 

Consider now an arbitrary permutation a in the per¬ 
manent expression perm A = n™=i which is 

actually a set {er(*)} containing the numbers from 1 to 
n in a certain order; further, be o' another permutation, 
obtained from er by changing the fc-th bit in all its compo¬ 
nents 0 ( 1 ) (written in their binary representation). Let’s 
analyse the effect of this bit flip. First, one should recall 
that, from its definition, 

( 8 ) 

with \a t j = 1 V i, j. Depending on the value of the 
fc-th bit of r^B, one has from 

U,b| fe = 1 => a iMi) = = 

= ~ '(i) = (i) (9) 

r i,B = 0 => ai jCT (i) = ( H{m)) ri (J ^ = 

= {H{m)) r . y({) = o ii<7 , (i) (10) 

If, as in the case of the hypotheses, an odd number of 
r^B has the fc-th equal to 1, in the product ri"=i a i,a'{i) 
an odd number of factors change their sign with respect 

t0 nr=i«^)> S ivin g : 

n n 

2=1 2=1 


( 11 ) 


5 


■+1 

+1 

+1 

+1 

+1 

+1 

+1 

+r 

+1 

-1 

+1 

-1 

+1 

-1 

+1 

-1 

+1 

+1 

-1 

-1 

+1 

+1 

-1 

-1 

+1 

+1 

-1 

-1 
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1 
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2 
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3 
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FIG. 1. Example of application of the criterion of Prop. [ 5 ] for a matrix A built of the rows {1, 2,3,3, 5, 6, 7,8} of H( 8). The 
table on the right summarise the application of the criterion: for each row i of the scattering matrix, which is the n line of 
7/(8), the binary expression is reported. The latter is actually the binary conversion of n — 1, since the binary enumeration 
of the rows must start from 0. The last line of the table reports the bitwise sum which, being not equal to zero in this case, 
indicates that perm A = 0. 


This means that for each permutation a there exist an¬ 
other one o', biunivocally associated to <7, for which © 
holds. Hence, in the sum over all the a of §f , half of the 
addends will have sign -1 and the other half + 1 , which 
implies perm A = 0. □ 

An example of application of this criterion is given in 
Fig. □ 

Corollary 5.1 Be A an m x m matrix, with m = 2 P > 
2, built by taking the rows {tT, r%,..., r TO } from H(m) 
(namely, the i-th row of A is the ri-th row of H(m) and 
rows may be repeated), and r\ t s © i~ 2 ,b © ■ ■ ■ © r mj B = 0. 
Build the matrix B such that all the rows are the same as 
those of A except the i-th, such i-th row being another ar¬ 
bitrary r\-th row of H(m), withr[ 7 ^ 7 + Then, B satisfies 
Prop. [d[ 

The condition t+b©^, 3 ®-. .©r mi B = 0 means that in 
the set of the binary representations {rps, r 2 ,s,..., r mj s} 
each bit recurs an even number of times with value 0 and 
an even number of times with value 1. The matrix B is 
built by removing from A its *-th row (which was the Ty-th 
of H(m)) and by replacing it with another arbitrary r'-th 
row of //(to). The binary representation K,b is different 
from r+B for at least one bit, say, the fc-th bit. Note that 
in the set {t+b, i' 2 ,b, ^m,s} such fc-th bit had the value 

0 for an even number of times, and the value 1 for an even 
number of times: now that we have changed 7 y with r[ 
the £:-th bit will have the value 0 for an odd number of 
times and the value 1 for (another) odd number times. 
This implies ri jB © r 2 ,B © ... © r' i B © ... © r n ,B 7 ^ 0 and 
the hypotheses for Prop. [5] are verified. □ 

Proposition 6 Be A an m x m, with m = 2 P > 2, built 
by taking the rows {tt, r 2 ,..., r n } from H(in) (namely, the 
i-th row of A is the ri~th row of H(m) and rows may be 
repeated). //r+Bffii+B©. • -©TVs =0, then perm A 7 ^ 0. 


Take the Laplace expansion of the permanent along an 
arbitrary i-th row: 


perm A = ajj ■ perm M t j 


( 12 ) 


where a 7 ;j is an element of A and Mij is the i,j minor 
of A. In other terms, (121 can be read as a dot product 


perm A = at ■ c 


(13) 


where a.j is the n-element vector given by the i-th row of 
A and c is the vector with elements Cj = perm M, j. 

Consider now m different matrices Ai , built by replac¬ 
ing the i-th row of A with the i-th row of //(to). For 
l = r.i one has A; = A r . = A, while all the other A; will 
differ from A by one row. When one calculates perm A; 
according to (12), the minors M i: j are always the same for 
every A/, because just the i-th row is changing. Hence, 
in (13) the vector c is always the same for every A;. The 


permanents of the different A/ can be interpreted as the 
projection of such c onto different vectors, given by the 
i-th row of //(to). 

It is important to note that c is a non-zero vector. The 
elements of this vector are permanents of matrices Mi j, 
which are squared (+ 1 ,- 1 ) matrices of order to — 1 and it 
has been shown [23] that if to = 2 P , then no matrices of 
order to — 1 exist with vanishing permanent. 

The rows of //(to) form a complete (orthogonal) ba¬ 
sis of R": a non-zero vector c has at least one non-zero 
projection on one of the vectors of the basis. We have 
already shown (Prop. 5.1) that, for A; defined as above 
with i ^ rt, one has perm A; = 0, i.e. the projection of 
con all the rows of //(to), except the r+th, is vanishing. 
It follows that the projection on the r.j-th row must be 
non-zero: this implies perm A ri = perm A 7 ^ 0. □ 
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8 

16 

32 

64 


n 

2 

| ~ 0.33 

4 

10 

= 0.4 

16 , 

36 ' 

- 0.44 

— ~ 0 47 

136 — 

256 ^ 

528 — 

0.48 

1,024 ^ 

2,080 — 

0.49 

0.5 
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24 , 
35 ‘ 

~ 0.69 

240 

330 

~ 0.73 

2 ’ 880 ~ n 74 

3,876 — ' 

39,168 ^ 
52,360 “ 

s 0.75 

574,464 ^ 
766,480 ^ 

: 0.75 

0.75 

8 




5,600 

6,435 

~ 0.870 

428,736 ^ q>-7a 

490,314 — U-^74 

53,829,888 

61,523,748 

~ 0.875 

9,309,189,120 

10,639,125,640 

~ 0.875 

0.875 


TABLE I. Fraction of suppressed states over the possible output states, when injecting a m-modes Sylvester interferometer 
with n photons in the first n inputs. The number of suppressed states has been calculated by checking the criterion of Prop. [7] 
for each possible output state. In the last column the estimation from the formula (141 is given for comparison. 
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6.434 

6.435 

> 0.999 

490,058 

490,314 

> 0.999 

61,458,212 

61,523,748 

~ 0.999 

10,622,348,424 

10,639,125,640 

~ 0.998 

0.998 


TABLE II. Fraction of suppressed states over the possible output states, when injecting a m -modes Sylvester interferometer 
with n fermions in the first n inputs. In the last column the asymptotic value 1 — ^ is given for comparison. 


Proposition 7 Consider a unitary transformation m = 
2 P modes U m = ^=7J(2 P ) with p = k + q, and an input 

state with n = 2 q bosons h = (1 + n ■ c,..., n + n ■ c), 
where 0 < c < ( 2 k — 1). The output state g = 
(gi, g 2 , ■ ■ ■ ,9n) is suppressed if and only if gi, B ® 92 ,b ® 

... ® g nj B 7 ^ 0 , being gi tB the binary representation of 
gi — 1 (i.e. the binary representation of the output mode 
number, starting the count from 0) truncated to the q least 
significant bits, and ® the bitwise sum (XOR operation). 

As a first thing we address the case m = n, with 
h = (1,..., n). Here the result comes directly from con¬ 
sidering that the probability of an output configuration 
g = (< 71 , g 2 , ■ ■ ■, g n ) is proportional to perm A where S is 
a matrix whose i-tli row is the g ^th row of H(n). Be¬ 
cause of Prop. [5] and [ 6 ] such permanent is vanishing if 
and only if g\^ B ® g 2 , B ® ... ® g n ,B 7 ^ 0 (where gys is the 
full binary expression of gi — 1 , composed of q bits) thus 
giving in this case the suppression of the corresponding 
output configuration. 

By exploiting Prop.[4j this result can now be extended 
to the case m = 2 q > 2 P = n and input states of the kind 
h = (1 + n ■ c,... ,n + n ■ c). In particular, the condition 
of Prop.[4]of considering the mode numbers modulo n, im¬ 
plies that a criterion of the kind gi, B ®g 2 ,B ®- ■ • ®gn,B 7 ^ 0 
can be applied if gi^ B is the binary representation of the 
mode index, truncated to the q least significant bits. □ 

To evaluate the fraction of output combinations that 
is suppressed we need to consider the set of all pos¬ 
sible output states {gi, g 2 ,..., g n } and estimate when 
their binary expressions {. 9 i,s, 92 ,b, ■ ■ ■, 9n,B }, truncated 
to the q least significant bits (for the arguments discussed 
above), satisfies G = gi,B®g 2 ,B®- ■ ■®9n,B 7 ^ 0. Adopt¬ 
ing an approach similar to that of Ref. [|7|, we assume 


that in such a set, a certain fc-th bit of the binary expres¬ 
sion g itB (consisting of q bits) can take the values 0 or 1 
with equal probability, independently from the values of 
the other bits. In other words, we assume that in each 
subset of states with a certain bit combination (for the 
bits other than the fc-th), the number of output states for 
which the fc-th bit is 0 is equal to the number of states 
for which that bit is 1 . 

Let’s now consider the possible values of the binary 
expression G, starting from its first bit. That bit is the 
result of the ® operation on n bits (the first bit of each 
gi, B ). If we consider the full set of possible outputs, such 
n bits will be 0 or 1 the same number of times. Thus, also 
the first bit of G will be 0 for half of the possible output 
states and 1 for the other half. Those states for which 
the first bit of G is 1 already satisfy G 7 ^ 0, so they are 
suppressed. For the other ones, they may be suppressed if 
the ® operation on other bits give 1. One then considers 
the second bit and with analogous arguments notes that 
it will be 0 for the half of the output states and 1 for 
the other half. One continues with the same procedure 

w bos 

up to the o-th bit. Hence, the overall fraction » r Bupp of 

states 

suppressed states will be \ (fraction of states which has 
the first bit of R equal to 1 ) summed to \ • | (fraction 
of states which has the first bit of R equal to 0 and the 
second equal to 1 ) summed to | ^ ■ \ (fraction of states 
which has the first and second bit equal to 0 and the 
third equal to 1) and so on. This gives: 

jybos Jl 1 pfc _ I ,, _ 1 
supp ^ \ - _£_ _ £_£ _ ™ 

^states ~ ^ 2 * 2 k n K 

Table H] reports the fraction of suppressed bosonic 
states for n < 8 and m < 64, compared with the re¬ 
sult of (fl4|). The latter expression approximates better 
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the actual value for increasing n or to. 


V. MULTIPLE FERMIONS 

The case of fermions is less significant, with respect to 
bosons, from a computational point of view; in fact the 
probability amplitude of the output configurations are 
proportional to the determinant of the scattering matrix 
(see Eq. ©)> which differently from the permanent can 
be calculated efficiently. However, investigating the sup¬ 
pression laws arising for this kind of particles enable a 
better understanding of the effects of statistics in multi¬ 
particle interference. 

Proposition 8 Consider a unitary transformation over 
to = 2 P modes U m = 2 P ) with p = k + 

q, and an input state with n = 2 q fermions h = 
(1 + n ■ c,... ,n + n ■ c), where 0 < c < (2 k — 1). The 
output state g = (gi, ■ ■ ■ ,g n ) is suppressed if and only if 
gi mod n = i mod n Mi £ [1, n]. 

Let’s consider, to begin, the case n = to, i.e. n 
fermions entering a n-mode interferometer one per each 
port. The only possible output state allowed from the 
Pauli principle is g = (1,... ,n), namely the state hav¬ 
ing one particle per mode, which can be written also 
as gi = i Vi £ [l,n]. Such condition is easily ex¬ 
tended to the case of a more generic input state h = 
(1 + n ■ c,... ,n + n ■ c) of n = 2 q fermions entering a 
to = n ■ 2 k interferometer through Prop. [4j becoming 
gi mod n = i mod n Mi e [l,n]. □ 

The number of allowed/suppressed output states can 
be evaluated considering that, in an interferometer with 
m = n ■ 2 k modes, the condition gi mod n = i mod n 
can be satisfied for 2 k different values of gi. Hence, the 
number of allowed output states is the number of se¬ 
quences of n numbers, each with 2 k possible values, i.e. 

2 k-n _ ^ m ^ n 

The possible output states of n particles on to modes 
are Abates = ( m+ ^ _1 ) = ( (m-’lyn! (combinations with 
repetitions). The fraction of allowed states, for large to, 
thus tends to[M]: 


(m\ n (to — 1)! n! 

( m \ 

. n n\ n! 

\nJ (to + n — 1)! 

\m-l) 

n n n n 


and the fraction of suppressed states is asymptotically 
equal to: 


K 


fer 


supp 


N si 


= 1 - 


N! 


fer 

allowed 


JV si 


1 - — (16) 
n n 


Table [TT] reports the fraction of suppressed fermionic 
states for n < 8 and to < 64, calculated over the all the 
possible n-particles states, together with the asymptotic 
estimation with (16). 




FIG. 2. a) Ratio between the average number of occu¬ 
pied modes for bosons and for distinguishable particles, for 
n = 2, 4, 8 particles entering an m-mode Sylvester interferom¬ 
eter, one per each of the first n modes, b) Probability dis¬ 
tribution of detecting an output state with a certain number 
of occupied modes, when entering with n = 8 indistinguish¬ 
able bosons (one per each input mode) in a m = 8 Sylvester 
interferometer (blue squares) or Fourier interferometer (red 
diamonds). Probability for distinguishable particles (black 
circles) in analogous interferometers is also reported. 


VI. DISCUSSION 


The predictions of the suppressions laws proved above 
can be compared with the general trends that usually 
distinguish the particles behaviour, depending on their 
statistics. While for fermions the compatibility of Prop.JH] 
with the Pauli principle is intrinsic in its same proof, the 
behaviour of bosons, which would be expected to show 
an enhanced tendency to bunch together, is more pe¬ 
culiar. Actually, similarly to Ref. ?], we observe that 
many bunching events are not enhanced but included in 
the suppression conditions of Prop. [7] The probability of 
full-bunching events (events with all n bosonic particles 
on the same output mode, over the to possible modes) can 
be calculated easily: the scattering matrix is composed 
of identical rows with half +1 f \frn and half — 1 /y/rn ele¬ 
ments; by inverting the sign of the columns with negative 
elements (operation which does not affect the permanent 
[25]), we obtain a matrix of all +1 /\/to, whose perma¬ 
nent equals n\/m m / 2 . Thus the probability of this event 

is = n -/ mm (f rom squaring Eq. 1) which is 

an enhancement of n\ with respect to the probability of 


















such event for distinguishable particles. This enhance¬ 
ment factor for full bunching events is indeed a general 
law for all unitary processes [3J. With regard to anti¬ 
bunching events (particles in all different ports), it is 
not difficult to observe that for n = m > 2 they are 
instead never suppressed for bosons: the scattering ma¬ 
trix for such an event would be a full Sylvester matrix, 
whose permanent is proved [2S3 to be non-vanishing by 
observing that for {r i, r^, ■ ■ ■, r n } = {1,2,..., n} one has 
ri,B © T 2 ,b © • • • © r n ,B = 0 and by exploiting Prop. [6] 
This also marks a difference with respect to Bell muh 
tiports, where such events are always suppressed for an 
even number of bosons 21 ' • 

An overall figure that may quantify the bunching be¬ 
haviour [7] is the average number M of occupied inodes 
at the output. The bosonic bunching tendency should 
reduce the number of occupied modes with respect to 
the case of distinguishable particles. Figure reports 
the ratio between the average number of occupied modes 
in the case of indistinguishable bosons and that in the 
case of distinguishable particles, for n = 2,4,8 parti¬ 
cles injected in the first modes of Sylvester interferome¬ 
ters with up to 32 modes. This ratio is always smaller 
than unity, confirming the bunching behaviour. How¬ 
ever, while from the previous discussion we know that 
the fraction of suppressed states is practically constant 
with increasing to, this ratio looks to approach one in the 
same limit: the larger number of available modes makes 
the particles more likely to exit on different ports, ap¬ 
proaching the classical probability on this aspect. Thus, 
the suppression law seems to be a stronger non-classical 
signature than the bunching behaviour itself, which is 
less evident with large to. 

The detailed probability of having a certain number 
of occupied modes is reported in Figure [2]o for the case 
n = to = 8, with the comparison of the distribution in the 
case of a Fourier interferometer. The three distributions 
are different in shape and it is evident the shift towards 
a smaller number of occupied modes for the two non- 
classical distributions. In this particular case the average 
number of occupied ports is ~ 4.1 for both Sylvester and 
Fourier interferometers in case of identical bosons and is 
~ 5.3 in case of distinguishable particles. 

As a further analysis, the asymptotic fraction of al¬ 
lowed fermionic states can be compared with the asymp¬ 
totic fraction of allowed bosonic states, to evaluate the 
strength of the suppression law in the two cases. For 
large to, one reads: 


Snowed _ ^ = n n-1 
^states n n n n 


1 1 

< 

n n 


j y b °s 
allowed 

N states 


(17) 


Interestingly, the suppression law seems to act more 
severely for fermions, for large to. Note that the calcu¬ 
lation in actually includes, for fermions, all the pos¬ 
sible multi-particle states, even those already forbidden 
from the Pauli principle itself. It can be observed, how¬ 
ever, that N^ tes = (™) x ( m+ ” _1 ) = Abates for large 
to, where N^l tes are all the possible fermionic states, i.e. 


states with all different output ports. Thus the inequality 
holds asymptotically also considering (for fermions) 
only the fraction of events that were not already sup¬ 
pressed by the simple application of the Pauli principle. 
It is worth reminding that, besides manifesting naturally 
for true fermions and bosons, the effects of the statis¬ 
tics can be simulated by proper entangled states 123111]. 
Thus, the laws here developed for the two kinds of par¬ 
ticles hold true for the corresponding entangled states. 

From a more applicative point of view, suppression 
laws such as the one presented here for Sylvester interfer¬ 
ometers may be exploited to test the indistinguishability 
of n-photon sources^- Further, they may be used, in 
the context of boson sampling experiments, to simulta¬ 
neously check the quality of the sources and of a possible 
reconfigurable device [18], that would perform the re¬ 
quired unitary that expresses the suppression law. The 
limitation of the laws for Sylvester interferometers to spe¬ 
cific values of n may seem, at first glance, quite disad¬ 
vantageous in a possible general case when it may be 
required to test a general n-photon source where n may 
not be a power of 2. However, one may envisage that 
the full n-photon interference could be used as an overall 
check of the source quality, but an accurate troubleshoot¬ 
ing of possible malfunctioning or imperfections requires 
different subsets of single photon sources to be tested sep¬ 
arately. To this purpose one may configure a device to 
implement a block-diagonal matrix, having for each block 
a different 2 p -mode (Sylvester) unitary to test separately 
different subsets of single-photons. 

An interesting perspective to test at the same time 
the indistinguishability of different couples of photons, 
without the need of reconfiguring the device is given by 
Prop. [3] which is characteristic of Sylvester matrices and 
does not hold for Fourier ones. In fact, in the two photon 
case, whichever couple of inputs is excited, an identical 
fraction of outputs is suppressed. This may be partic¬ 
ularly useful to check single-photon sources for scatter¬ 
shot boson sampling [22(: there, several heralded single¬ 
photon sources are used, each coupled to a different in¬ 
put port, and multi-photon states (with one photon per 
port, but multiple ports excited) are generated randomly. 
Couples of photons on random input ports are generated 
efficiently in such a setup and by comparing the detected 
output events with the predictions of Prop. [3] the indis¬ 
tinguishability of all the possible couples of sources may 
be conveniently tested. 


VII. CONCLUSION 

In conclusion, we have proved a necessary and suffi¬ 
cient criterion for the suppression of many output com¬ 
binations when n = 2 q particles are injected in certain 
inputs of a linear interferometer with to = 2 P modes im¬ 
plementing a Sylvester matrix. While both the bosonic 
and fermionic cases have been studied, the result is par¬ 
ticularly significant for bosons, whose output distribution 
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is hard to compute in the general case. Therefore, this 
suppression law may be exploited for the use of Sylvester 
multiports as benchmark devices for the indistinguisha- 
bility of multiple single-photon sources or the assessment 
of the overall quality of reconfigurable interferometers. 

This study has also shown that comprehensive laws 
that describe the output multi-photon distribution of 
multi-port interferometers on the basis of the symme¬ 
try of the implemented matrix are not limited to Fourier 
ones. Indeed, further investigations could pursue the def¬ 
inition of similar criteria for wider class of matrices, thus 


giving greater insight on the features of multi-particle 
interference. 
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